Pseudogap behaviour in Bi2Ca2SrCu208: Results of Generalized Dynamical 

Mean-Field Approach 
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Pseudogap phenomena are observed for normal underdoped phase of different high-Tc cuprates. 
Among others Bi2Sr2CaCu208-5 (Bi2212) compound is one of the most studied experimentally. 
To describe pseudogap regime in Bi2212 we employ novel generalized ab initio LDA+DMFT+Ek 
hybrid scheme. This scheme based on the strategy of one of the most powerfuU computational 
tool for real correlated materials: local density approximation (LDA) + dynamical mean-field the- 
ory (DMFT). Here conventional LDA-I-DMFT equations are supplied by an additional (momentum 
dependent) self-energy Ek in the spirit of our recently proposed DMFT-(-Ek approach, account- 
ing for pseudogap fluctuations. In the present model Ek describes non-local correlations induced 
by short-ranged collective Heisenberg-like antiferromagnetic spin fluctuations. The effective single 
impurity problem of the DMFT is solved by numerical renormalization group (NRG). Material spe- 
cific model parameters for effective x'^ — orbital of Cu-3d shell of Bi2212 compound, e.g. the 
values of intra- and interlayer hopping integrals between different Cu sites, local Coulomb interac- 
tion U and pseudogap potential A were obtained within LDA and LDA+DMFT. Here we report 
theoretical LDA+DMFT-f Ek quasiparticle bands dispersion and damping, Fermi surface renormal- 
ization, momentum anisotropy of (quasi) static scattering, densities of states, spectral densities and 
angular resolved photoemission (ARPES) spectra accounting for pseudogap and bilayer splitting 
effects for normal (slightly) underdoped Bi2212 (5=0.15). We show that LDA-|-DMFT-|-Ek suc- 
cessfully describes strong (pseudogap) scattering close to Brillouin zone boundaries. Our calculated 
LDA-|-DMFT-|-Ek Fermi surfaces and ARPES spectra in presence of the pseudogap fluctuations are 
almost insensitive to the bilayer splitting strength. However, our LDA-calculated value of bilayer 
splitting is found to be rather small to describe experimentally observed peak-dip-hump structure. 
Results obtained are in good semiquantitative agreement with various recent ARPES experiments. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.27+a, 71.30.-(-h, 74.72.-h 



I. INTRODUCTION 

Pseudogap state is the major anomaly of the normal 
state of copper oxides, commonly believed to be most 
relevant for the understanding of the physical nature of 
high-Tc superconductivitjsi. 

During the last decade experimental techniques of an- 
gular resolved photoemission spectroscopy (ARPES) has 
made a large progress. Sate-of-the-art high-Tc test com- 
pound for ARPES is Bi2Sr2CaCu208-5 (Bi2212) system. 
Thus a lot of experimental ARPES data is available for 
Bi2212 (for reviews seeS). Several major experimental 
characteristics are derived from ARPES data like, for 
example, Fermi surfaces (FS), quasiparticle band disper- 
sions and damping, and even self-energy lineshapes^. A 
number of interesting physical anomalies were discovered 
in normal underdoped phase of Bi2212: pseudogap for- 
mation, shadow bands and bilayer splitting of FS^ . These 
phenomena abounds in theories and there arc still no def- 
inite point of view about their physical origin. It is be- 
lieved that all of these are quite relevant for the physics 
of high-temperature superconductors. The problem is 
much complicated by strong electronic correlations ever 
present in these compounds, making the standard band 
theory and Fermi liquid approaches doubtful. 

In this work we show that taking into account short- 
range antiferromagnetic fluctuations resulting in pseu- 



dogap formation together with bilayer splitting effects 
is enough, in principle, to describe abovementioned 
experiments. To this end we employ novel hybrid 
ab initio LDA-|-DMFT-|-Ek computational scheme^i^i^. 
From one side this scheme inherits all the advantages 
of LDA-|-DMFT6''^.i^iaiio, i.e. the merger of first prin- 
ciple one-electron density functional theory within local 
densty approximation (DFT/LDA) -*^^'^^ and the dynam- 
ical mean-field theory (DMFT) for strongly correlated 
electron o-'^^d^d^d^d'^ . On the other side our scheme al- 
lows one to account for non-local correlation effects intro- 
ducing momentum dependent external self-energy pre- 
serving conventional DMFT equations^i^. To solve the 
effective single impurity problem of the DMFT we em- 
ploy here the reliable numerical renormalization group 
approach (NRG)^^'^^. 

Such combined scheme is particularly suitable to de- 
scribe electronic properties of real high-Tc materials at 
finite doping in the normal state. First, all material 
specific model parameters for physically relevant effec- 
tive — orbital of Cu-3d shell are obtained from 
LDA computations. Second, undoped cuprates are an- 
tiferromagnetic Mott insulators with U ^ W {U — 
value of local Coulomb interaction, W — bandwidth 
of non-interacting band), so that correlation effects are 
very important. Thus at finite doping (up to optimal 
doping) cuprates are typical strongly correlated metals. 
To this end DMFT stage in our computational scheme 
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takes these strong electronic correlations into account. 
To adopt LDA+DMFT to study the "antiferromagnetic" 
scenario of pseudogap formation in cupratea^ i^"'^^'^^ k- 
dependent self-energy Ek describing non-local correla- 
tions induced by (quasi) static short-ranged collective 
Heisenberg-like antiferromagnetic (AFM) spin fluctua- 
tions is included^iiS^. 

Recently we applied DMFT+Sk approach to inves- 
tigate formation of pseudogap for strongly correlated 
metallic regime of single-band Hubbard model on the 
square lattice'^'^''^. At present there are several inde- 
pendent methods aimed to describe non-local effects be- 
yond standard DMFT. Similar results about pseudo- 
gap formation in the 2d Hubbard model were already 
obtained within two-particle self-consistent approach^, 
cluster DMFT extensions, such as the dynamical 
cluster approximation (DCA)^^ and cellular DMFT 
(CDMFT)2^^^^, CPTia^i^ and via the model of 
two interacting Hubbard sites self-consistently embed- 
ded in a batb^. The EDMFT was employed to demon- 
strate pseudogap formation in the DOS due to dynamic 
Coulomb correlations'^'^. Important progress was also 
made with weak coupling approaches for the Hubbard 
model^ and functional renormalization group)^"^. In 
several papers pseudogap formation was described in the 
framework of the t-J model^l. A more general scheme 
for the inclusion of non-local corrections was also formu- 
lated within the so called GW extension to DMFT^^i^. 
Dynamical vertex approximation to study Mott-Hubbard 
transition in presence of non-local antiferromagnetic cor- 
relations was proposed^. Chain-DMFT extension was 
used to investigate breakup of the Fermi surface near 
Mott transition for quasi Id Hubbard model''^. 

The paper is organized as follows. In section |TT] 
we present a short introduction into an ab initio 
self-consistent generalized combined LDA-I-DMFT+Ek 
scheme to account for short-rang AFM correlations. Sec- 
tion mil contains Bi2212 material specific information: 
LDA calculated band structure and details on some 
model parameters calculations. Results and a discus- 
sion of LDA-I-DMFT+Sk computations for Bi2212 are 
presented in the sections IIVI and |Vl 



II. COMPUTATIONAL METHOD 

A. Introduction of the length scale into DMFT: 
DMFT+Sk approach 

To introduce spatial length scale into conventional 
DMFT methodi^ii^'i^ii^'i^ recently we proposed gener- 
alized DMFT-|-Ek approacb^i^. The major assumption 
of our approach is that the lattice and Matsubara "time" 
Fourier transform of the single-particle Green function 
can be written as 



G(a;,k) 
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where S(ci;) is the local self-energy of DMFT, while 
I]k(w) is some momentum dependent part. Interfer- 
ence effects between these parts are neglected. Advan- 
tage of our generalized DMFT-t-Sk approach is addi- 
tive form of self-energy in Eq. ([T])'^'^'^. It allows one 
to keep the set of self-consistent equations of standart 
DMFT-i^iiiii^ii^iil. However there are two distinctions. 
First, on each DMFT iteration we recalculate corre- 
sponding k-dependent self-energy Ek(/Lt, w, [E(w)]) within 
some (approximate) scheme, taking into account interac- 
tions with collective modes or order parameter fluctu- 
ations. Second, the local Green's function of effective 
impurity problem is defined as 



M ^ 



N ^ iuj + fi- e(k) - S(w) - Ek(w) ■ 



(2) 



at each step of the standard DMFT procedure. 

Eventually, we get the desired Green function in the 
form of ID), where E(a;) and Sk(w) are those appearing 
at the end of our iteration procedure. 

To calculate Sk(^) for an electron moving in the 
random field of pseudogap fluctuations (assumed to be 
(quasi) static and Gaussian, which is valid at high 
enough temperature o^^'^^ ) with dominant scattering mo- 
mentum transfers of the order of characteristic vector 
Q = (7r/a,7r/a) (a - lattice spacing) of AFM fluctuations 
( "hot spots" model^ ) , we use the following recursion pro- 
cedure proposed in Refs i^^i^^''*^ 



where 



S„(c^,k) = 



I]k(cc;) =S„=i(cc;,k), 



s{n) 



(3) 



iuj + ^ - S](w) - e„(k) + invnK - E„+i(w, k) 

(4) 

The quantity A characterizes the pseudogap energy scale 
and K — is the inverse correlation length of short 
range SDW fluctuations, e„(k) = £(k+Q) and u„ — 
I'^^k+ql + bk+ql o*^^ while e„(k) = e(k) and 
Vn = l^^kl + l^kl ^'^^ even n with w^^^(p) determined 
by usual momentum derivatives of the "bare" dispersion 
£(k), while s(n) represents a combinatorial factor, deter- 
mining the number of Feynman diagrams^^'^^. 

For the (Heisenberg) spin structure of interaction 
with spin fluctuations in "nearly antiferromagnetic 
Fermi-liquid" (spin-fermion (SF) model of Ref.^^) spin- 
conserving scattering processes obey commensurate com- 
binatorics, while spin-flip scattering is described by dia- 
grams of incommensurate type ( "charged" random field 
in terms of Refill). In this model combinatorial factor 
s{n) acquires the following formal 



sin) 



^ for odd n 
J for even n. 



(5) 



+ Ai-e(k)-S(tj)-Sk(t^)' 



(1) 



Obviously, with this procedure we introduce an impor- 
tant length scale ^ not present in conventional DMFT. 
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Physically this scale mimics the effect of short-range 
(SDW) fluctuations within fermionic "bath" surround- 
ing the effective Anderson impurity of the DMFT. We 
expect that such a length-scale dependence will lead to a 
kind of competition between local and non-local physics. 

Though we prefer to consider both parameters A 
and ^ as phenomenological (to be determined by fitting 
experiments)*, one can in principle calculate these from 
microscopic model at hand. For example, using the two- 
particle self-consistent approach of Refsi^i^ with the ap- 
proximations introduced in Refsj ^^'^^ , we derived within 
the standard Hubbard model the following microscopic 
expression for A* 

A2^C/2^^^^1?^<(n.T-naf >, (6) 

where we consider only scattering by antiferromagnetic 
spin fluctuations. The different local quantities - to- 
tal density n, local densities n^i, and double occu- 
pancy < Uiiriii > - can easily be calculated within the 
standard DMFT—. A detailed derivation of ([6]) is pre- 
sented in the Appendix B of Ref.**. Corresponding micro- 
scopic expressions for the correlation length ^ = k^^ can 
also be derived within the two-particle self-consistent 
approach^M^. However, we expect these results for 
to be less reliable, because this approach is valid only 
for relatively small (or medium) values of U/t, as well 
as for purely two-dimensional case, neglecting quasi-two- 
dimensional effects, obviously important for cuprates. 
Actually, our calculation experience shows that all the 
results obtained below are rather weakly dependent on 
the values of ^ from the experimentally relevanii interval 
of (5 H- 10)a. 

B. Bilayer splitting effects: LDA+DMFT+Ek 
formulation 

To perform ab initio calculations for Bi2212 system we 
employ LDA-^DMFT strategy proposed in Refs-^^i^i^ii^. 
Necessary bare band dispersion for effective physically 
relevant Cu-3d — orbital in a tight-binding repre- 
sentation is 

e(k) = —2t {cos k^a + cos kya) (7) 
—At' cos k^a cos kyQ 
—2t" (cos 2kxa + cos 2kya) 
—2t"' (coskxacos2kya + cos2kyacoskya). 

Here t, t', t" , t'" are hopping integrals within first four 
coordination spheres. Tight-binding equation for inter- 
layer dispersion is taken in the form 

t i^{k) — ^ [cos krj-a — cos kyaY (8) 

given in Ref.*"* with bilayer splitting equal to 2tj_. 

Since account of bilayer splitting (BS) effects in Bi2212 
requires essentially two-band model we introduce bare 



Hamiltonian in reciprocal space as the following matrix 
over (bonding and antibonding) band indices 

The local Green's function is now also a matrix 

G(^) = ^ E 0^ - H(k) - (E(^) + S]k(cc'))i) . (10) 

where we have assumed self-energies to be diagonal. In 
the following we would like to keep the DMFT part of the 
problem just a single-band task. This can be achieved 
taking the diagonal element of (fTO|) 

^_ 1 G-i(u;,k) 
GM-]^l.(G-i(c.,k))2-(i^(k))2' (11^ 

where G{uj, k) is given by ([T]). This local Green's function 
G{oj) (which includes additive self-energy contributions) 
determines now our effective single Anderson impurity 
problem. One should remark since we work with the 
single-band problem there is no need for double counting 
correction between LDA and DMFT.— 



III. LDA BAND STRUCTURE OF BI2212 AND 
EFFECTIVE MODEL PARAMETERS 

The Bi2212 compound has tetragonal bcc crystal lat- 
tice with symmetry space group 14/mmn*^'^'^. Main 
structural motif for this compound is two Cu02 layers 
displaced close to each other in the unit cell. Using crys- 
tal structure data of Ref.*^ we performed LDA calcula- 
tions of electronic band structure within the linearized 
muffin-tin orbital (LMTO) basis set^. Obtained band 
structure is in agreement with the one of Ref.*^. 

In Fig. [1] one-electron LDA band dispersion along BZ 
symmetry lines for Bi2212 is shown. Gray lines corre- 
spond to all-band Hamiltonian. To extract physically in- 
teresting partially filled — orbital of Cu-3d shell 
Wannier functions projecting method^ in the LMTO 
framework^*^ was applied. Corresponding dispersion of 
effective — y^ orbital is displayed in Fig. [1] as a black 
line. 

To set up the LDA-f-DMFT+Sk lattice problem ^ 
one needs to calculate transfer integrals t, t' , t" , t'" and 
t for tight-binding expressions ([7]) and ([H]) . On the basis 
of Wannier function projecting^ we performed computa- 
tion of corresponding hopping integrals with its LMTO 
realization^-. Obtained values for intra- and interlayer 
hybridization between — y^ orbital of different Cu-sites 
are listed in the Table I. Values of t, i', t", t'" we present 
are somewhat larger than those extracted from ARPES 
experiments^. On the other hand our value of t^ is much 
smaller than experimental one t^''=0.083 eV— . At the 
same time our calculated value of t_L is in good agreement 
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with other band structure results reported^. Taking into 
account large difference between t± and further we 
provide LDA+DMFT+Sk results for both these values. 

The value of local Coulomb interaction U for x'^ — 
orbital was obtained via constrained LDA method^. To 
screen this x"^ — orbital we used the rest of the Cu-3d 
shell of our selected site, neighbouring inplane Cu sites 
and also Cu sites from closest Cu02 layer. The value 
found is U=1.51 eV (Table I). 

Pseudogap potential A (see Eq. was obtained as 
described in Ref4. using LDA+DMFT(NRG) to calculate 
set of occupancies entering (|6]) (instead of DMFT(QMC) 
used in Ref.^). For given values of hopping integrals and 
U value with hole doping level S = 0.15 our A equals 
0.21 eV. The value of correlation length ^ is always taken 
to be equal to 5 lattice constants which is a typical exper- 
imental value^. Temperature comes through NRG part 
of our scheme and is always taken to be '^255 K. This 
completes the set of necessary model parameters to start 
LDA+DMFT+Ek computations for Bi2212 (see SecHIl). 



IV. RESULTS AND DISCUSSION 
A. Bi2212 LDA+DMFT+Ek densities of states 

Density of states (DOS) is calculated as 

N{lj) = --lmG{uj) (12) 

TT 

where G{uj) is defined by Eq. [TT] analytically continued 
to real frequencies. In Fig. [2] we display LDA+DMFT 
and LDA+DMFT+Sk DOS for effective x^ - orbital 
of Cu-3d. It is clearly seen that pseudogap fluctuations 
lead to formation of the pseudogap in DOS within 0.2 eV 
from the Fermi level. In our model this pseudogap is not 
tied to the Fermi level and it is not very pronounced for 
parameter values used here for Bi2212. It is also easy to 
find out that for all our DOS curves BS effects are most 
pronounced on the top of the van Hove singularity, which 
is about around -0.2 eV below the Fermi level (see also 
inset of Fig.^for details). Namely we calculate DOS for 
LDA value of BS 0.03 eV (gray curve) and experimental 
BS value of 0.083 eV (black hue). For the latter case BS 
effects are obviously stronger. Dashed curves correspond 
to LDA-fDMFT results for two different values of bilayer 
splitting. For LDA+DMFT+Sk DOS (solid curves) it is 
observed that BS effects become less pronounced (but 
still can be seen for the case of ^^^=0.083 eV). This is 
caused by the decrease of the life-time due to pseudogap 
fluctuations. Also van Hove singularity becomes slightly 
narrower here due to self-energy effects. Note that the 
shape of the pseudo gap in the DOS almost does not 
depend on BS effects. 



B. Bi2212 LDA+DMFT+Sk qusiparticle 
dispersions and damping 

For the case of flnite temperature and interaction val- 
ues we define quasiparticle dispersions via maxima posi- 
tions of corresponding spectral functions A{uj, k) 

A{uj,k) = --lmG{uj,k), (13) 

TT 

where G{uj,]i) is defined by an expression under the sum 
in (jlip . analytically continued to real frequencies, with 
self-energies and chemical potential ^ calculated self- 
consistently as described in Sec. Ill Al 

In Figs. [3] and [1 we present LDA+DMFT and 
LDA-|-DMFT-fSk quasiparticle bands dispersions 
(crosses) for Bi2212 effective x'^ — orbital of Cu-3d 
shell along the symmetry lines in the Brillouin zone (BZ) 
for t± and t']^^. Background shows quasiparticle damp- 
ing given by the imaginary part of additive I](w) -|-Sk(w) 
self-energy. The more intensive shade corresponds to 
the larger damping. In case of standard LDA+DMFT 
computations, neglecting non-local corrections (Fig. [3]), 
one can clearly see that the damping is uniform over all 
BZ. This is due to local nature of conventional DMFT. 
Quasiparticles are well defined in narrow light region 
around zero energy (Fermi level). 

When we introduce spatial inhomogeneity into DMFT 
bath within the LDA+DMFT+Sk approach the damp- 
ing appears to be much stronger and consequently non- 
uniform as seen in Fig. S) Again quasiparticles are well 
defined close to the Fermi level. But now the contour plot 
of Im[I](w) + Ek(w)] self-energy (damping) clearly shows 
so-called "shadow band" which looks like the quasiparti- 
cle band mirrored around the zero energy. In Fig. [T] we 
can also see pseudogap formation around X point. In our 
case shadow band is formed due to short-ranged AFM 
fiuctuations. Close to X point BS effects are most pro- 
nounced. One can see that maxima of A{uj, k) belonging 
to the "shadow band" region are conserved only rather 
close to the X point, further away these maxima vanish 
due to large damping. In the middle of MG direction we 
observe preformation of AFM insulating gap in the cross 
point of quasiparticle and and "shadow" bands. 



C. Bi2212 LDA+DMFT+Ek spectral functions 

To plot spectral functions A{u!,'k.) (fT5|) we choose k- 
points along the 1 / 8-th part of the "bare" Fermi surface 
within the first quadrant of the Brillouin zone for given 
lattice spectra and filling. In Fig. [5] corresponding spec- 
tral functions for different strength of bilayer splitting are 
shown. 

Close to the to nodal point (upper curve) spectral func- 
tion in Fig. [5] has the typical Fermi-liquid behaviour, 
consisting of a rather sharp peak close to the Fermi level. 
Going to the antinodal point (lower curve) fiuctuations 
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becomes stronger and shift the sharp peak out of the 
Fermi level down in energy. Simultaneously with the 
growth of fluctuation strength damping also grows, so the 
peak becomes less intensive and more broad. In the vicin- 
ity of the "hot-spot" (black line) the shape of A{uj,'k) 
is completely modified. Now A{uj,'k) becomes double- 
peaked and non-Fermi-liquid~like. Directly at the "hot 
spot", A{ui,'k) has two peaks (second one is much less 
intensive) situated symmetrically around the Fermi level 
and splitted from each other by ^ 1.5A^^'^^. 

For the case of (right panel of Fig. [5]) behaviour 
is similar to one for tj_ (upper panel of Fig. [5]). How- 
ever now bilayer splitting strength is big enough to be 
resolved. So the peak-dip-hump structure^ is formed on 
the edges of pseudogap. 



higher binding energies as we move from nodal to antin- 
odal region. Displacements of theoretical and experimen- 
tal peaks on the left panel of Fig. [7] are in good quanti- 
tative agreement. However theoretical peaks are always 
a little bit sharper and narrower. Take notice that the 
left panel demonstrates no BS effects. On the right panel 
of Fig. [7] we found slightly better agreement of intensi- 
ties due to bigger BS value t^j^^. But for the k- values 
between the "hot spot" and the antinodal point there 
we have some lack of spectral weight close to the Fermi 
level. Also for these k-values we observe some reminis- 
cence of the bilayer splitting. After all one can infer that 
BS effects do not change the line shape of our ARPES 
spectra significantly and for both cases we obtain rather 
satisfactory agreement with the experiment. 



D. Bi2212 LDA+DMFT+Ek ARPES spectra 

Knowing A(uj, k) (jl3p we are now of course in a posi- 
tion to calculate angle resolved photoemission (ARPES) 
spectra, which are the most direct experimental way to 
observe pseudogap in real compounds. For that pur- 
pose, we only need to multiply our results for the spec- 
tral functions with Fermi function at temperature 255 K. 
The resulting LDA+DMFT+Sk ARPES spectra are pre- 
sented in Fig. [6l Again these spectra are drawn along 
1/8 of non-interacting FS from antinodal (lower curve) 
to nodal point (upper curve). At the antinodal point 
we find well defined (sharp) quasiparticle peak close to 
the Fermi level. Moving towards the antinodal point 
the damping (widening) of this quasiparticle peak and 
its shift to higher binding energies are observed. Such 
behaviour is typically obtained experimentally^. To de- 
scribe peak-dip-hump splitting resolved in experiment^ 
we take i^''=0.083 eV^^. Indeed for i^^'' we get pro- 
nounced peak-dip-hump structure similar to experimen- 
tal one^. It is recognized that our LDA-calculated t± 
is several times smaller and can not provide adequate 
description of the peak-dip-hump structure for ARPES 
data. Notice that the intensity of antibonding branch 
is higher than for bonding one. It is opposite in the 
experiment. We attribute this difference to the matrix 
elements effects which are not taken into account in the 
present work. 

In Fig. El we show comparison of LDA-|-DMFT-|-Sk 
ARPES spectra and experimental one from ReffS^ for 
Bi2212 measured along the Fermi surface. Here spec- 
tral functions displayed in Fig. [5] are multiplied with the 
Fermi function at experimental temperature T=140 K 
and convoluted with Gaussian to simulate experimental 
resolution of 16 meV,^ All theoretical ARPES curves 
after multiplication and broadening are normalized to 1. 
Left and right panels of Fig. [7| correspond to the theo- 
retical data for t± and t^^^ values. Both figures demon- 
strate semiquantitative agreement of our theoretical re- 
sults with the experiment. Common trend for both pan- 
els is the damping of quasiparticle peak and its retreat to 



E. Bi2212 LDA+DMFT+Ek Fermi surface 

In the following we characterize renormalized Fermi 
surfaces (FS) by intensity plots of spectral density at zero 
frequency A{u! — 0, k) (which for the free-electron case 
just follow the "bare" Fermi surface). 

In the Figs. [5] and [5] we display thus defined 
LDA-hDMFT and LDA-KDMFT-^Sk Fermi surfaces for 
Bi2212. LDA-hDMFT FS has the LDA shape, as it 
should be within DMFT (see Fig. [5]). Slight broaden- 
ing close to the borders of BZ is because of BS effects. 
Non-zero width of FS (in contrast to LDA) comes from 
finite damping due to interaction and temperature. For 
LDA+DMFT+Sk FS (see Fig. O one can see signifi- 
cant "destruction" effects in the vicinity of the antinodal 
point induced by pseudogap fluctuations. From compar- 
ison of upper and lower panels of Fig. [9] one can conclude 
that for strongly correlated case BS effects alone are not 
enough to describe experimentally observed FS "destruc- 
tion" close to the borders of BZ and formation of "Fermi 
arcs" around the nodal point, as observed in ARPES 
experiments^. We found that FS shape is rather insen- 
sitive to BS strength since pseudogap fluctuations are 
much stronger than bilayer splitting and hide it. Though 
BS at the BZ boundaries slightly amplifies pseudogap ef- 
fects. Thus the account of pseudogap (AFM) fluctuations 
seems to be necessary to describe experimental picture. 



F. Bi2212 LDA+DMFT+Ek anisotropy of static 
scattering 

Strong anisotropy of (quasi) static scattering was 
observed in Bi2212 system in ARPES experiments in 
Refs.^i i^^i^^ and attributed to scattering by planar 
impurities^i^. Here we show that this effect can be 
naturally explained by (quasi) static scattering by pseu- 
dogap fluctuations. 

Our LDA+DMFT+Ek calculated (quasi) static scat- 
tering defined as a(k) = E(0) + Sk(0) is plotted in 
Fig. [ini together with experimental data of Refsj ^*^'^^ . 
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Here k-points are taken along 1/8 of non-interacting FS. 
We detect our results to mediate experimental data of 
Refs.— while the difference between latter remains it- 
self unexplained. 

In our opinion, anisotropy of (quasi) static scatter- 
ing a(k) naturally follows from anisotropic renormaliza- 
tion of electronic spectrum due to pseudogap fluctua- 
tions, which directly follows from our "hot spot" like 
model^ii^. 

Despite overall behaviour is analogous to one obtained 
in the experiment there is a need for further studies of 
possible relevance of matrix elements effects in ARPES, 
as well as that of additional scattering by random static 
impurities^. 



V. CONCLUSION 

Present investigation is aimed to describe pseudogap 
regime of high-Tc cuprate Bi2Sr2CaCu208-5 (Bi2212) 
from first principles. For this purpose we employ novel 
generalized ab initio LDA-I-DMFT-I-Ek hybrid scheme. 
This scheme based on the strategy of most powerfull com- 
putational tool for real correlated materials: local den- 
sity approximation (LDA) -I- dynamical mean-field the- 
ory (DMFT). Here we supply conventional LDA+DMFT 
equations with an additional (momentum dependent) 
self-energy Ek in the spirit of our recently proposed 
DMFT-I-Sk approach. "External" self-energy Sk is cho- 
sen here to describe non-local dynamical correlations in- 
duced by short-ranged collective Heisenberg-like antifcr- 
romagnetic spin fluctuations (in static Gaussian approxi- 
mation of RefS i^^'^^ ). Necessary Bi2212 material specific 
model parameters for effective x'^ — orbital of Cu-3d 
shell, e.g. the values of intra- and interlayer hopping inte- 
grals, local Coulomb interaction U and pseudogap poten- 
tial A were calculated within LDA and LDA-^DMFT. On 
the basis of LDA-t-DMFT+Ek computations we obtain 
densities of states, spectral functions A(uj, k) which allow 
one to visualize quasiparticle bands dispersion and damp- 
ing, Fermi surface (FS), anisotropy of static scattering 
a(k) and ARPES spectra accounting for pseudogap and 
bilayer splitting effects for normal (slightly) underdoped 
Bi2212 ((5=0.15). It is found that on the DOS level BS 
and pseudogap effects are separated in energy and hardly 
affect each other. We showed that LDA+DMFT+Sk de- 
scribes strong scattering at Brillouin zone boundaries as 
pure manybody effect. LDA-fDMFT-l-Sk Fermi surface 
in presence of the pseudogap fluctuations is almost in- 
sensitive to the BS strength. Thus the BS effects alone 
are not enough to describe the Fermi surface destruction 
(though amplifies it) and additional source of electron 
scattering is required (for example, AFM short-range 
fluctuations). The only place where BS effects play sig- 
nificant role is formation of the experimentally observed 
peak-dip-hump structure in ARPES spectra. To this end 
the LDA-calculated value of bilayer splitting is found to 
be rather small to describe this effect. Results obtained 



are in good semiqualitative agreement with various re- 
cent ARPES experiments. 

At present there are several alternative points of view 
on the possible explanation of Fermi surface destruction, 
formation of shadow Fermi bands etc. Recently the anal- 
ysis of the effect of three-dimensionality on the ARPES 
spectra was presented for Bi2212 in Refi^. It was shown 
that in a quasi-2D system, the weak /c^-dispersion can 
lead to Fermi surface maps similar to those observed in 
the experiment. This FS broadening mechanism does not 
have the manybody origin. The authors of Refi^ have 
shown that shadow Fermi surface in Bi2212 can be in- 
terpreted as an intrinsic feature of the initial electronic 
spectrum arising from bulk, orthorhombic distortions lo- 
cated primarily in the BiO planes, but most definitely 
felt throughout the three-dimensional crystal. All these 
effects are not considered here thus remain for further 
investigations. Apparently, in a real system these mech- 
anisms combine with those described above leading to a 
complete picture of electronic structure of Bi2212. 
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Figures 




FIG. 1: Calculated within DFT/LDA Bi2212 band disper- 
sions (gray lines) and effective x'^ — band of Cu-3d shell 
obtained by projection on Wannier functions (black lines). 
The Fermi level corresponds to zero. 
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FIG. 2: LDA+DMFT (dashed lines) and LDA+DMFT+Ek 
(solid lines) densities of states for Bi2212 for LDA-calculated 
value of ix=0.03 eV (black) and experimental value of 
^ea:p_g Qgg ^gj.g^y') (Coulomb interaction U=1.51 eV, filling 
n=0.85, pseudogap potential A=0.21 eV, correlation length 
^ = 5a). Inset shows magnified region around the Fermi level. 




FIG. 3: LDA+DMFT quasiparticle bands for Bi2212 (crosses) 
along BZ high symmetry directions for LDA-calculated value 
of tx=0.03 eV (upper panel) and experimental value of 
^ea:p_g Qgg ^jQ.^gj. panel) (Coulomb interaction U=1.51 
eV, filling n=0.85). Zero of background (which is -1 / ■n\xaS{ijj) 
- local DMFT self-energy) corresponds to zero damping. 
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FIG. 4: LDA+DMFT+Ek quasiparticle bands for Bi2212 
(crosses) along BZ high symmetry directions for LDA- 
calculated value of t_L=0.03 eV (upper panel) and experimen- 
tal value of t'j^^=0.083 eV (lower panel) (Coulomb interaction 
U=1.51 eV, filling n=0.85, pseudogap potential A=0.21 eV, 
correlation length ^ = 5a). Zero of background (which is 
-l/7rIm[E(tj) + Ek(tj)] - additive local and "pseudogap" self- 
energies) corresponds to zero damping. 
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FIG. 5: LDA-I-DMFT-I-Ek spectral densities for Bi2212 along 
of noninteracting FS in 1/8 of BZ for LDA-calcuIated value 
of t±=0.03 eV (upper panel) and experimental value of 
^ea:p_Q Qgg ^^^.^gj- panel) (Coulomb interaction U=1.51 
eV, filling n=0.85, pseudogap potential A=0.21 eV, correla- 
tion length ^ = 5a). 
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FIG. 6: LDA+DMFT+Ek ARPES spectra for Bi2212 along 
of noninteracting FS in 1/8 of BZ for LDA-calculated value 
of t±=Om eV (left) and experimental value of fj'^ =0.083 
eV (right). (Coulomb interaction U=1.51 eV, filling n=0.85, 
pseudogap potential A=0.21 eV, correlation length ^ = 5a). 
Corresponding spectral function A(lu, k) is multiplied with 
Fermi function at T~255K (the temperature of NRG calcu- 
lations) . 
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FIG. 8: LDA+DMFT Fermi surfaces for Bi2212 within 1/4 
of BZ {kx,ky in units of w/a) for LDA-calculated value 
of tx=0.03 eV (upper panel) and experimental value of 
^e3:p_Q Qgg ^jQ.^gj- pajiei) (Coulomb interaction U=1.51 

eV, filling n=0.85). 
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FIG. 7: Comparison of LDA-|-DMFT-|-Ek ARPES spectra 
(solid lines) for Bi2212 along of noninteracting FS in 1/8 of 
BZ for LDA-calculated value of tx=0.03 eV (left) and ex- 
perimental value of t1^^=0.083 eV (right) with experimental 
ARPES (Ref.^4) (circles). (Coulomb interaction U=1.51 eV, 
filling n=0.85, pseudogap potential A=0.21 eV, correlation 
length ^ = 5a). Corresponding spectral function A{ijj,'k) is 
multiplied by Fermi function at T=140K (the temperature of 
experiment) and broadened with Gaussian to simulate exper- 
imental resolution of 16 meV (Ref.^''). 
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FIG. 9: LDA+DMFT+Ek Fermi surfaces for Bi2212 within 
1/4 of BZ {kx,ky in inits of n/a) for LDA-calculated value 
of tx=0.03 eV (upper panel) and experimental value of 
^ea;p_Q ggg ^^.^ygj. panel) (Coulomb interaction U=1.51 
eV, filling n=0.85, pseudogap potential A=0.21 eV, correla- 
tion length ^ — 5a). 
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FIG. 10: Comparison of experimental and theoretical 
LDA+DMFT+Ek static scattering a(k) = -l/7rImE(0) + 
Ek(0) for Bi2212 within 1/8 of BZ for LDA-calculated 
value of tx=0.03 eV (light curve) and experimental value of 
^e3:p_g Qgg ^jjark curve) (Coulomb interaction U=1.51 eV, 
filling n=0.85, pseudogap potential A=0.21 eV, correlation 
length ^ = 5a). 
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TABLE I: Calculated energetic model parameters for Bi2212 
(eV). First four Cu-Cu inplain hopping integrals t, t' , t" , t'" , 
interplain hopping value t± , local Coulomb interaction U and 
pseudogap potential A. 
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t' 


t" 


t'" 


t± 
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-0.627 


0.133 


0.061 


-0.015 


0.03 


1.51 


0.21 



